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of interspecies interaction. Here, we include intraspecific competition to 
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regime are some examples of these characteristics. From the stability of the steady 
state solutions, three phases are obtained: (i) extinction of one species (Cause 
transition), (ii) their coexistence and (iii) a forbidden region. Full analytical 
solutions have been obtained for the considered ecological regimes. The time 
transient allows one to defined time scales for the system evolution, which can 
be relevant for the study of tumor growth by theoretical or computer simulation 
models. 
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1. Introduction 

The ecological community is composed by a complex interaction network from 
which the removal of a single species may cause dramatic changes throughout the 
system. The interactions between species only became better formalized in population 
dynamics with the Lotka-Volterra equations, in the 1920's pQ. This simple theoretical 
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model explains prey-predator behavior and others range of natural phenomena, for 
instance, the oscillatory behavior in a chemical concentration in chemistry or in fish 
catches in ecology, cells of the immune system and the viral load (in immunology) [3j 
econophysics [3] etc. 

Besides the predation interaction described by Lotka-Volterra equations, there 
are many other different kinds of interaction taking place between biological species. 
If species unfavor each other, for instance, when two species occupy the same 
ecological niche and use the same resources, there is competition. If species favor 
each other, as in pollination/seed dispersion by insects, there is mutualism - also 
called symbiosis. If only one of these species is independent of the other, there are 
two possibilities: amensalism, if the considered species has a negative effect on the 
other; and commensalism, otherwise. As an example of amensalism, consider that an 
organism exudes a chemical compound as part of its normal metabolism to survive, 
but this compound is detrimental to another organism. An example of commensalism 
is the remoras that eat leftover food from the shark. Finally, when species do not 
interact at all there is the so-called neutralism. 

The Verhulst-like two-species population model incorporates limit environmental 
resources, logistic growth of one species in the absence of the other, and interspecific 
interaction. It is described by the following equations [US]: dN±/dt — k 1 A^ 1 (I — 
Nx/Ki + a 1 N 2 /K 1 ) and dN 2 /dt = n 2 N 2 (l - N 2 /K 2 + a 2 N 1 /K 2 ), where Ni > 0, 
Ki and Ki > are the number of individuals (size), net reproductive rate, and the 
carrying capacity of species i (= 1,2), respectively. The term —KiNf/Ki represents 
the competition between individuals of the same species (intraspecific competition), 
and — K1N1N2J K\ represents the interaction between individuals of different species 
(interspecific interaction) . The carrying capacity represents the fact that both species 
are resource-limited. In fact, K\ represents the restriction on resources that comes 
from any kind of external factors, but the ones relative to species 2, similarly for 
K 2 . To use non-dimensional quantities in the Verhulst-like two-species model, write 
Pi = Ni/Ki > 0, for i = 1,2. Time is measured with respect to the net reproductive 
rate of species 1, r = Kit > 0. Here, we restrict ourselves to the case n t > 0. The scaled 
time is positive since we take the initial condition as to = 0. Moreover, the two net 
reproductive rates form a single parameter p — k 2 /ki > 0, fixing a second time scale to 
the system: r' = pr = n 2 t. The parameters a\ and a 2 represent the interaction degree 
and also the kind of interaction between species. The non-dimensional population 
interaction parameters are given by t\ = a.\K 2 jK\ and e 2 = a 2 K\j K 2 , which are not 
restricted and represent the different ecological interactions. With these quantities, 
the Verhulst-like two-species model become: 

^7 = Pi[l-Pi +£1^2] = f(pi,P2) (1) 

= -P2 + £2Pi] = g(pi,P2) ■ (2) 

Contrary to p, which has no major relevance to this model (since we consider only 
p > 0), the product 6162 plays an important role, so that e\e 2 < means predation; 
€162 = means commensalism, amensalism, or neutralism; and eie 2 > means either 
mutualism or competition as depicted in the diagram of Figure [l] 

The novelty of the present work is two-fold. We give a new interpretation of the 
Verhulst-like two-species model, since we do not restrict the interaction parameter 
€i, as it is usually done in previous studies [51 O [7]- This absence of restriction 
on the interaction parameters allows us to display many different ecological regimes 
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(competition, predation, and mutualism) in the same phase diagram, which are 
obtained by the steady state solutions [HI M EH EH E3 EH- Nevertheless, we point 
out the existence of a survival/extinction transition, as well as a transition for a 
forbidden region in the mutualism regime. Also, based on the result of Ref. |14) . 
we have been able to find the complete analytical solution for this model. From 
the analytical solutions, the transient time allows one to establish characteristic 
times [HI US H71 [HI HSJ EDI EU [22] which can be relevant, for instance for tumor 
growth models [531 [Ml US]. Also transient solutions may be used as a guide to 
validate numerical simulation algorithms, giving to Monte Carlo time steps a real time 
scale [26 . Furthermore the analytical solutions may be used in stochastic theoretical 
modeling of two-species systems to give insight when averages are done [27] ■ 

The text is organized as follows. In Sec. [2j the steady state solution of the 
interspecific competition in the Lotka-Volterra equations is presented. These solutions 
correspond to the stable ecological regimes in the parameter space diagram, where we 
point out the existence of the survival/extinction transition and of a non-physical 
region in the mutualism regime. We show that although our model is simpler, it is 
able to reproduce the four scenarios found in more complete tumor growth models. 
In Sec. [3j the analytical solutions for the trivial case neutralism and the non-trivial 
ones amensalism and comensalism are presented. Next, the full analytical solutions 
are obtained for the mutualism, predation, and competition regimes. Our conclusions 
are described in Sec. [H 



2. Steady state solutions 

To obtain the steady state solution p\ — pi(r —¥ oo) and p 2 = P2^J oo) 
of Equation ([!]) and Equation p]), we have to impose dpi/dr = dpifdr = 0, 
which implies /(p*,^) = QipliP^) = an( i leads to p*(l — p* + eiP 2 ) = and 
pZO- ~~ P*2 + e 2Pi) = 0. One has the following trivial ("t"), semi-trivial ("st"), and 
non-trivial ("nt") pairs of solutions: (i) p\ t = and p 2 t — 0; (ii) p\ st — 1 and 
P*2,st = or Ptst = and P*2,st = 1 and (iii) p* lnt = (1 + ei)/(l - eie 2 ) and 
P%. nt = (1 + e 2)/(l — £i£2)- These asymptotic solutions characterize the system 
according to their stability. The stability matrix, also called community matrix £Q, 
{ d f d f \ 

is: A{p\,p%) = I J 2 I . The steady state solutions p* and p* 2 are 

stable if the trace and the determinant of the community matrix are negative and 
positive, respectively. One has: Tr[A(pJ,p2)] = 1 + P + {p<^2 - 2)p* + (ei - 2p)p* 2 and 
Det [A(p\ , p* 2 )] = p{ 1 + p\ [e a - 2(1 + e 2 pl )] + p* 2 [ei - 2(1 + e^)] + Ap\p* 2 } . One must 
analyze the possible cases given p\ and p* 2 ■ 



2.1. Stability analysis 

Let us start with the stability analysis of the trivial solutions of p* t = and p 2 t — 0, 
which means extinction of both species (synnecrosis). One has: Tr[A(0,0)] = 1 + p 
and Det[A(0,Q)] = p. Since p > , DetL4(0,0)] > but Tr[A(0,0)] > 1. The pair 
of trivial solution is not stable anywhere in the parameter space, so synnecrosis never 
occurs in our model. 

The semi-trivial solutions are given by p\ st = 1 and p 2 st — or p* st — and 
P2 st — 1 and they mean that one of the species is extinguished. Considering the 
species 1 extinction, one has: TrL4(0, 1)] = 1 + t\ — p and DetL4(0, 1)] = — p(l + ei). 
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For these solutions to be stable, it is necessary that e\ < — 1, regardless of the p value. 
A similar analysis leads us to conclude that species 2 extinction is stable only for 
£2 < -1. 

The non-trivial solutions p\ <nt = (l + ei)/(l — e^) andp2.„ t = (1 + £ 2)/(l — £i£2) 
lead to: T^W 1M ,p* 2M )] = [1 + e 1 + (1 + e 2 )p}(e 1 e 2 - 1) and Det[A(pI |Tll> p5 )Trt )] = 
— (1 + ei)(l + e 2 )p/{e\e 2 — !)• On one hand, if £162 < 1, the denominator is positive 
and the numerator of p\. n t and p2,nt must vanish or be positive. From the condition 
Pint ^ 0) this solution is only stable if e\ > e\ = — 1; otherwise, p* lt = is the 
stable solution. This produces a transition from the regime where species 1 coexists 
with species 2 to the regime where species 1 is extinguished (Cause transition). The 
same transition occurs for the parameter e 2 . O n the other hand, if e\t 2 > 1: the 
denominator is positive and the numerator of p% >n t an d P2,nt must be non-negative. 

(c) 

From the condition, pi nt > 0, this solution is only stable if t\ < e\ = —1; otherwise, 
p\ t — is the stable solution. This produces a transition from the regime where 
species 1 coexists with species 2 to the regime where species 1 is extinguished. The 
same transition occurs for e 2 - 

From the stability criteria, species can coexist only if ei > —1 and e 2 > — 1. 
According to the values of e\ and e 2 , the various ecological regimes may present a 
stable non-trivial solution, as in Figure [T] 

2.2. Phase diagram 

In Figure [TJ the stable steady state solutions of Equation ([!]) and Equation |2]) are 
represented in the parameter space. This diagram presents the reflexion symmetry 
about e 2 = t\ and summarizes our findings: the coexistence phase and one species 
extinction phase can be seen. These phases span on different ecological regimes, which 
means that different ecological interactions may lead to the same phase. 

For ei > — f, e 2 > — f, the non-trivial solutions are stable. They span on four 
considered ecological regions, around the origin of the phase diagram [see the circle in 
Figure]!]. For mutualism, the first quadrant of the diagram of Figure[lJ as £162 - > 1 — , 
the mutual cooperation conducts to unbounded growth of both populations, so that 
Pi,nt ~ (1 — ei£2)~^ diverges with the exponent /3 = 1 (see Figure [2]). The region 
e 2 > 1/ei is forbidden; since p* nt < 0, it does not have ecological reality. 

The Figure [2] also tell us that there is a transition between extintion and 
coexistence regimes. Keeping e 2 fixed, this transition occurs at t\ = ef = — 1. A 
Taylor expansion of the non-trivial solution allow us to write p* nt = (ei — £i)/(l + 

e 2 ) + o((ei — e1) 2 ^j , whith means that near (the critical point) species 1 linearly 

goes extinct, that is p* nt ~ (ei — ef). The critical exponent related to the order 
parameter is /? = 1 Analogous process happens to the species 2. 

For ei < — 1 or e 2 < — 1 and £162 < 1, the remaining regions of the diagram are 
characterized by the stability of the semi-trivial solutions. For e\ < — 1 and £2 < — 1, 
differently from all the other regions in the phase diagram, the steady state solutions 
depend on the initial condition. There is a separatrix for the initial conditions, in this 
region pQ. 
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Figure 1. Diagram of ecological interaction between two species according to the 
pair (ei,e2). In this diagram, each quadrant represents one kind of interaction: 
I:mutualism; II and IV: predation ; III: competition. The abscissa and ordinate 
represent either ammensalism or commensalism. The origin represents neutralism. 
The non-trivial solutions (ei > —1, €2 > —1, and e\£2 < 1) correspond to the 
coexistence phase. In this phase mutalism, comensalism, amensalism, predation, 
competition, and neutralism can occur, for instance, around and inside the dashed 
circe. The complementary region is characterized by an extinction phase and a 
forbidden region. The extinction phase (ei < —1 and/or t2 < —1) reveals a region 
for ei < and €2 < where, contrary to the other cases, the steady state solutions 
depend on the initial condition. For ei > 1 and t2 > and £2 > 1/ei, there exists 
a forbidden region with non biological reality (negative number of individuals). 



2.3. Cancer modeling 

In cancer therapy, tumor-selective replicating viruses offer remarkable advantages 
over conventional therapies and are a promising new approach for human cancer 
treatment. An oncolytic virus is a virus that preferentially infects and lyses cancer 
cells. Theoretical models of the interaction between an oncolytic virus and tumor 
cells are attainable using adaptations of techniques employed previously for modeling 
other types of virus-cell interaction 0. A Lotka-Volterra like model that describes 
interaction between two types of tumor cells (the cells that are infected by the 
virus and the cells that are not infected but are susceptible to the virus so far as 
they have the cancer phenotype) and the immune system has been presented by 




Figure 2. Plot of steady steady solutions p* t = (1 + ei)/(l — ei€2) and 
P2 „ t = (1 + £2)/(l — ei«2) of the model [Equation |TJ and Equation Q] as a 
function of e\ , keeping €2 = 1 fixed. In these plots one sees two phase transitions. 
The first one is the extinction-coexistence (predation) transition, which occurs at 
ei = —1. Near £j, the critical point, the species 1 go to extinction by a linear 
form, that is p* nt ~ (ei — ej). The second one is the transition from coexistence 
(mutualism) to a forbidden region (with no ecological reality), which one with 
critical exponent (9 = 1. 



Wodarz 23, 24] , This model can be written as [25]: dx/dr = x [1 — (x + y)] — f3xy; 
dy/dr = 72/ [1 — (x + y)]—/3xy—5y, where, x and y are the population of uninfected and 
infected cells respectively, 7 is the ratio between infected and uninfected cells growth 
rates, j3 is related to the interaction parameter between uninfected and infected cells 
and (5 is related to the rate of infected cell killing by the virus. Using Equation ([!]), 
Equation ^ and making: t\ — —(1 + 0), e 2 = {fih — 1) and 7 ^> S, one retrieves 
Wodarz model. Besides being simpler, our model presents all the qualitative behavior 
shown by Wodarz model, i.e. four different scenarios for the asymptotic states: (i) 
absence of infected cells, (ii) absence of uninfected cells, (hi) coexistence of both types 
of cells and (iv) dependence on initial conditions (all cells infected or uninfected). 
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3. Model analytical solution 

Here we present the analytical solution of the Verhulst-like Lotka-Volterra model. We 
start presenting the known solutions when one of the interacting parameter vanishes. 
Next we present the solution for the complete model. 



3.1. Vanishing of one interaction parameter (t\ei = 0) 

This section is restricted to the particular case £162 = 0, where one or both interaction 
parameters vanish. This corresponds to the axis of the parameter space C2 x e\ 
[see Figure IT]. Thus, three ecological regimes are allowed in this specific situation: 
amensalism: t\ = and 62 < (species 2 extinction, if £2 < — 1 and species coexistence 
otherwise) or e 2 = and e\ < (species 1 extinction, if ex < —1 and species coexistence 
otherwise); neutralism: e± = e 2 = 0; and comensalism: e± > and £2 = or e\ = 
and e 2 > 0. 

In these cases, one can obtain a simple full analytical solutions of Equation (JlJ 
and Equation |2]). Below we address each case in more detail. 

3.1.1. Neutralism Consider a special case where each population grows independent 
from the other. This ecological regime is represented by Equation (JlJ and Equation 
([2]), with ei — e 2 = 0, leading to two independent Verhulst models: dpi/dr = pi[l— pi] 
and dp2/(pdr) — p2[l — p?]. The solutions, with different time scales and parameters, 

for each species are [H]: Pi( T ) = l/[l+(Pi~o — l) e ~ T ] andp2(T") = l/[l + (P2^o — l) e_pr ]' 
where pi t o — pi(0) is the initial condition for species i = 1,2. 

The Verhulst solutions are driven by different characteristic times r = Kit and 
t' = pr = n 2 t, respectively. For r ^> 1, so that pr ^> 1, the asymptotic behaviors 
Pi = Pi(°°) = 1 an d P*2 = P2{oo) = 1 are obtained, so that species end exploring all 
the available environmental resources. The case 62 = in Figure [3] shows the dynamics 
of the population given by the Verhulst solutions. For k 2 > «i; i.e., p > 1, species 2 
grows more rapidly than species 1, given the same initial condition. For K2 < n±; i.e., 
p < 1, the inverse occurs. 



3.1.2. Comensalism and amensalism Consider that two species interact asymmet- 
rically. For instance, consider that individuals of species 1 are unaffected by species 
2, although, individuals of species 2 are adversely affected by species 1. This is the 
amensalism regime. The commensalism regime has the same structure as amensalism, 
except that one species is favorably affected by the other. These interactions can be 
mathematically represented by the following equations: 
dpi(r) 



dr 
dp 2 {r) 



: Pi {t)[1- Pi {t)] (3) 



= P2(t)[1-P2(t) + c 2Pi (t)] , (4) 
par 

where e 2 is negative for amensalism and positive for commensalism. 

In this kind of interaction, species 1, described by Equation j3j|, follows the 
Verhulst model, whose solution is P\{t) = 1/[1 + {pi \ — l)e~ r ] . The dynamics 
of species 2 follows the time-dependent Verhulst-Schaefer model [Equation Q], whose 
solution is [T4"] : 

r^r = r\ + 1 ( 5 ) 
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p 1 (r) = z / dr' Pl (r') = \ _ , ^ = ln[l + Pl , (e T - 1)] . (6) 



where the mean relative size of species 1 up to t is 

i r r dy 

la i + (pr,o-!)e 

Using the Verhulst solution for p\ (t) and put the Equation ^ in Equation (J5j) , one 
obtains: 

1 _ e- T [l + {-l + eJ)p lfi ] 2 



Pa(r) Pi,o{l + [-l + e r (l + e 2 )]p^ } 

e-P T [! + (-! + e r ) Pi,o]~ 62PT (1 + e 2 Pi,o - P2,o) 
£>2,0 + e 2 Pi : oP2,o 

The plots of P2(r) for several e 2 values are depicted in Figure [3j 



(7) 




Figure 3. Plots of the solution of Equation given by Equation |7j, for 
different values of interaction parameter E2 and p = 1. The comensalism regime 
is obtained for 62 > 0, where the asymptotic value, due to the other species, 
is greater than unity. The neutralism regime is retrived for t2 =0. The 
amensalism regime is obtained for £2 < 0, where the asymptotic value does not 
vanish (meaning species coexistence) for £2 > — 1 and vanishes (meaning species 
extinction) for €2 < — 1. 

The steady state solutions of Equation ^ and Equation Q are, respectively: 
pi = pi(oo) = 1 and p* 2 = ^2(00) = 1 + e 2 , if e 2 > —1 or vanish otherwise. One sees 

(c) 

that e 2 = — 1 is a critical value that separates two distinct phases: e 2 < — 1, where 
species 2 is extinguished; and e 2 > —1, where species 2 coexists with species 1. The 
former case occurs in the amensalism regime, while in the latter one it may occur in 
the amensalism (e 2 < 0), neutralism (e 2 =0), or in the comensalism (e 2 > 0) regimes. 
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The same conclusions are valid for £2 = and t\ ^ 0. One finds the same 
aviors am 
extinguished. 



behaviors and a critical point ef^ = — f, so that similarly for e\ < species 1 is 



3.2. Mutualism, competition and predation (ex £2 ^0) 

In the following, we deal with the case £162 ^ 0, which addresses mutualism, 
competition, and predation. If 

• £162 > 0, each species has the same kind of influence on the other. This 
corresponds to either the competition or the mutualism regime. The following 
regimes occur: 

— £1 > and £2 > 0, mutualism, which corresponds to the first quadrant of the 
space parameter phase space, restrict to the region £ 2 < 1/ej (see Figure [lj 

— £1 < and £ 2 < 0, competition, which corresponds to the third quadrant of 
the space parameter phase space (see Figure [I]) . 

• If £i£2 < 0, the predation regime occurs, which belongs to the second and fourth 
quadrants of the parameter space (see Figure [T]) . 

(c) 

For £2 > and e\ < 0, there is species coexistence for £1 > e\ ' = — 1 and species 1 

extinction for ei < ej'. 

These ecological regimes are special cases of Equation ([I]) and Equation ([2]), whose 

solutions can be worked out to have the form: 

1 1 , e- T [ 1+ei ^M](l-Pi,o + eiP2,o) 

— (°) 



Pl(r) 1 + £iP 2 (t) Plfl + £lPl,0P2,0 

where P2(t) is given by Equation ([5]), and the relative populations sizes mean values 
up to instant r are Pi(t) — fl cIt'pi(t')/t and p 2 ( r ) = Jo dr' P2(t') / t . 

Using Equation ([5]) in Equation ([8]), we obtain a quadratic equation for Pi(t), 
eliminating its dependence on P2 (r). In fact, we can write Pi{t) as dependent only 
on the initial condition and l^ 7 ")- The population size P2(j) behaves analogously. 
Thus, the coupling between the two population sizes is given only by the mean values. 
The solutions Equation ^ and Equation ^ are presented in Figure [4] for the three 
regimes where £i£2 7^ 0. As r — > 00, the steady state solutions of Equation ([T]) and 
Equation ^ are reached. 

With the analytical solution, one can access the transient behavior of a given 
population. Transient dynamics can be an important aspect of the coexistence of 
predators and preys, and also of competitors [TS]. Studies of outbreaks (insects or 
diseases) focuses greatly on the transient dynamics [THl H7J [T51 [THl HOI 1211 122- In 
tuberculosis treatment for example, it can reveal important aspects beyond asymptotic 
states such as how drug resistance emerges [28]. In the present study Figure [i] 
illustrates the importance of the transient in the time evolution of the species densities. 
Consider the case of competition, for r < 5 species 1 population is greater than species 
2, however the steady state solution for this system is just the opposite (species 2 
population is greater than species 1). In this case, a simple steady state analysis 
would not be coherent with reality if the observation time scale is not appropriate. In 
other words, if the observed system has not yet reach equilibrium, the steady state 
analysis can be misleading. 

Notice that considering e\ = in Equation ^ and Equation ([5]), one retrieves 
the Verhulst solution and Equation (|5j), and Equation which correspond to the 
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amensalism, neutralism, and comensalism regimes. In this way, these evolution 
equations can be seen as a general solution that is valid for all kinds of interaction 
regime. 




Competition 



S, 



-0.5 



C,=-0.1 



Mutualism 



S,=0.5 
E,=0.] 



Predation 
S,=0.5 
E 2 =-0.5 



Figure 4. Plots of the evolution equations Equation (|8]l and Equation j5), 
which are solutions of Equation ([TJ and Equation pi. The initial condition is 
Pl,0 = P2,o = 1/100 and p = 1. The species interaction parameters ei£2 ^ so 
that for competition: ei = —1/2 and t2 = —1/10; mutualism: ei = —1/2 and 
€2 = —1/10; and predation: ei = 1/2 and €2 = —1/2. 



4. Conclusion 

The simple model we addressed here illustrates that one can interpret the interaction of 
a two species system at several levels. From the interaction parameters ei and €2, which 
act at the individual level of species, one is able to tell about the different ecological 
regimes, classified in a higher level according to the product of the interaction 
parameter £162. If it vanishes, one or two species are independent from each other. 
If £162 > 0, one has either mutualism (both positive) or competition (both negative). 
For ei£2 < 0, one has predation. A collective level is obtained from the stability 
of the steady state solution, from where one obtains three phases: extinction of one 
species (e^ < —1) (synnecrosis is not a stable phase in our model), species coexistence, 
and a forbidden phase (e2 > 1/ei)- Although the studied model has been considered 
in several isolated instances, our study reveals the very general aspect of a simple 
mathematical set of equations, which represents very rich ecological scenarios that 
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can be described analytically. In this manuscript we focused on a Vcrhulst term 
for the population growth in the Lotka-Volterra equation. All the results presented 
here can be extended using a more general growth model, for instance the Richards' 
model [221 1301 EI] and implications of such generalizations will be detailed in a brief 
future. 
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